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Abstract 

This paper shows that, a quasi-localization of wavefunctions in toroidal compactifications 
with magnetic fluxes can lead to a strong suppression for relevant Yukawa couplings, and 
it is applicable to obtain tiny neutrino masses. Although it is known that magnetic fluxes 
lead to a Gaussian profile of zero-modes on a torus and that can yield a suppressed coupling 
in higher-dimensional supersymmetric Yang-Mills (SYM) theories, the largest (diagonal) 
entry of Yukawa matrices is always of 0(1). In this paper, we propose a way to induce an 
absolutely tiny global factor of Yukawa matrices. In two SYM theories defined in different 
dimensional spacetime, their bifundamental representations must be localized as a point 
in some directions. Overlaps of such point-like localized wavefunctions and Gaussian zero- 
modes give a global factor of Yukawa matrices, and it can be a strong suppression factor or 
a usual 0(1) factor, corresponding to their distance. Our numerical analysis shows that it is 
possible to obtain a suppression strong enough to realize the tiny neutrino masses without 
a hard fine-tuning. Furthermore, we propose a concrete model of the tiny neutrino Yukawa 
couplings in a magnetized SYM system. 
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1 Introduction 


Extra dimensional space is an attractive candidate for new physics beyond the standard model 
(SM). This has been being studied very actively since Kaluza-Klein theory was proposed to 
unify the forces of nature. Although our circumstances has got more complicated than those 
at that time, such a unified theory is expected to contain the extra dimensional space which 
is an origin of many flavors and the forces of the SM. Indeed, superstring theories which are 
candidates for the unified theory is defined in ten-dimensional spacetime, and the existence 
of extra dimensional space is indispensable for realizing the flavor structure in those theories. 
Higher-dimensional supersymmetric Yang-Mills (SYM) theories are expected as well-motivated 
higher-dimensional field theories since they appear in low-energy limits of superstring theories. 
They are also attractive from the phenomenological viewpoint even if we do not consider super¬ 
string theories. The reason is that they are the most simple theories containing supersymmetry 
(SUSY) which is another great candidate for the new physics as well as the extra dimensional 
space, and the SM can be embedded into the SYM theories with a large gauge group. 

Toroidal compactihcations with magnetic fluxes are able to yield the three forces and the 
flavor structure of the SM in superstring theories and higher-dimensional SYM theories[TJ [2]. 
Magnetic fluxes on two-dimensional (2D) tori lead to generations of the quarks and the leptons, 
and the Higgs fields as degenerate zero-modes. Zero-mode equations for these fields can be 
solved analytically thanks to simplicity of tori, and zero-mode wavefunctions are obtained as 
solutions of the equations. We see from the wavefunctions that the zero-modes are quasi- 
localized on the magnetized torus with a Gaussian profile. When zero-modes are localized far 
away from each other, their four-dimensional (4D) effective coupling is suppressed[3]. Thus, 
magnetic fluxes on tori can induce hierarchical Yukawa couplings, and the hierarchical masses 
and mixing angles of the quarks and the leptons could be potentially obtained on the magnetized 
torus. Indeed, a semi-realistic pattern of them was realized in Ref. [4]. 

The tiny neutrino mass is a remaining mystery of the SM. Although the magnetized toroidal 
compactihcations can yield suppressed Yukawa couplings, the largest diagonal entry of Yukawa 
matrices is always estimated to be of 0(1), and another mechanism to realize the tiny masses 
is needed. An elegant solution for explaining the tiny masses is given by introducing heavy 
Majorana mass terms for the right-handed neutrinos, which is a so-called seesaw mechanism[5]. 
However, there is a problem in toroidal magnetized compactihcation of SYM theories. Such 
a Majorana mass term cannot be obtained straightforwardly in the theories. We are forced 
to assume that such mass terms might be generated by nonperturbative effects[6] [7] and/or 
higher-dimensional operators. Although there is a wide variety of the seesaw mechanisms, they 
all require the presence of additional mass terms. Thus, this paper aim to propose a simple way 
to realize tiny Yukawa couplings for the neutrinos in the magnetized toroidal compactihcations. 

This paper is organized as follow. Section 2 is devoted to review the magnetized toroidal 
compactihcations in SYM theories. Zero-mode wavefunctions on the magnetized torus and 
their leading Yukawa couplings are shown. We give a main idea and result of this paper in 
Section 3. A numerical analysis is also performed to verify consistency with the tiny neutrino 
masses. In Section 4, we propose a concrete model with a specihc configuration of magnetic 
huxes on the tori in a SYM system, where our mechanism to yield the tiny neutrino Yukawa 
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couplings will work correctly. We conclude with a discussion of prospects in Section 5. 

2 Toroidal compactification with magnetic fluxes 

This section gives an overview of magnetized toroidal compactihcations on the basis of six¬ 
dimensional (6D) SYM theories compactihed on a two-dimensional (2D) torus, following Ref. [6] 
(See also | 8 J for a review of magnetized toroidal compactihcations). In this section, we use 
complex coordinates z and a complex 2D vector held A z defined as 

- = ^ {x 1 + rx 2 ) , z = {z)* , 

A z = — - -(r * At — A 5 ), A z = (A z )^ , 

imr 

where (a: 1 , a: 2 ) are the real torus coordinates and (A 4 , A 5 ) are extra dimensional components of 
a (real) 6D vector held. Periodicity of this torus is expressed by identifications of z ~ z + 1 
and z ~ z + r. Complex structure of this torus is denoted by r. This is dehned with the torus 
radius R in the torus metric as 


ds 2 = g a bdx a dx b = 2 h dzdz , 


where 

j=(27ri?) 2 ( R ) T h — 2 (2nR ) 2 . 

On this torus, we consider the following gauge potential containing magnetic huxes and 
Wilson lines, 

A z = -- (Mz + 

lmr 

The magnetic huxes and Wilson lines are represented by (N x iV)-matrices in U(N ) SYM 
theories. Since this paper concentrates on Abelian forms of magnetic huxes and Wilson lines, the 
matrices take nonvanishing values only for their diagonal entries. Note that the nonvanishing 
entries of magnetic huxes are forced to be integer because of the Dirac’s quantization condition. 

We discuss gauge symmetry breaking due to these magnetic huxes and Wilson lines. In 
the case of U( 2) SYM theory, this U( 2) gauge symmetry is broken down to 17(1) x 17(1) when 
the two diagonal elements of magnetic huxes (Wilson lines) are chosen to be different values. 
This is extended to the case of 17(77) SYM theories, and the gauge symmetry is broken as 
17(77) —» Ua U ( N a), where U(N a ) is an unbroken gauge subgroup of 17(77). 

We consider zero-mode equations for a two-component spinor held on this magnetized torus. 
Fields living in extra dimensional space are decomposed into the zero-mode and multiple 
Kaluza-Klein modes, and we focus on the zero-mode in the following. The two-component 
spinor is denoted by 
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and then, their zero-mode equations are given by 


d z ^ + + [A g , 4> + \ = 0, 

= 0 . 


In the assumption of gauge symmetry breaking U(N ) —>■ U(N a ), a bifundamental repre¬ 
sentation (N a , Nb) of two remaining gauge subgroups U(N a ) x U(Nb ) contained in tp± is denoted 
by tJj±. Their zero-mode equations are given by 


d z + 
d z - 


7 r 

21m r 

7T 

21m r 


Z + C a b) 

{M a bZ + Cafe) 


V- 


o, 

o, 


(1) 

( 2 ) 


where M ab = M a — M b and Cab = Ca~ Cb, and M a and Ca are relevant entries of M and C- 

These equations were studied in detail in Ref. [2]. Either the zero-mode equations for ip^ 
or gives a well-dehned wavefunction, which is related to the sign of magnetic fluxes M ab . In 
the case of M ab > 0, a bifundamental has normalizable solutions in the zero-mode equation 
m, and then, zero-modes of ip 0 // are eliminated. There appear only the zero-modes of 'ip 0 ^ with 
a degeneracy given by M ab . In contrast, the zero-mode equation ([2D yields the zero-modes of 
'ip ab when M ab < 0, and its degeneracy is given by |M a b|. This is a kind of chirality projection 
to yield a 4D chiral spectrum like the SM. Furthermore, we can identify these degenerate zero- 
modes with the generations of the SM particles. \M ab \ = 3 unit of magnetic flux leads to the 
three generations of quarks and leptons. 

Their 4D effective coupling constants are of our main interest because they would be Yukawa 
couplings of the SM which determine the masses and mixing angles of the quarks and leptons. 
The coupling constant of 4D effective held theories is given as an overlap integral of wavefunc- 
tions on the torus. Thus, we are required to study the zero-mode wavefunctions obtained as 
solutions of Eqs. ([ID and (f2j) . The zero-mode wavefunctions are obtained analytically, which are 
expressed by using Jacobi-theta functions. We use an index I ab to label the \M ab \ degenerate 
zero-modes, that is, I ab = 1, 2,..., \M ab \. The wavefunction of the I ab -th zero-mode r ip I p b for 
M ab > 0 is then given by 



( z ,t) 


0 Iab,Mab (Z + Cab/M ab , T ) , 


J\J"m gitiM zlm z /Im r 




'I/M 

0 


(Mz, Mr) , 


(3) 


where the Jacobi-theta function is defined by 


■& 


a 

b 


(Y t) = J2 


giri(a+£) 2 T ^2iri(a+£) (u+b) 


Normalization factors are determined as 


(4) 



et ge 1 ^ = 


5 


Ui 
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Figure 1: The squared absolute values of wavefunctions of the three zero-modes |'0+ 6 | 2 given 
in Eq. (J3]) with ( M ab = 3, ( ab = 0, r = i) are shown, where we set Rez = 0 and the horizontal 
axis corresponds to the other direction of torus. 


and we find 


JVm = 


21m t\M\ 

A 2 


1/4 


(5) 


where A denotes the area of torus. For M ab < 0, the wavefunction of ip T J lb is given by 

(z + C ab / Mab ,r)Y . 


These zero-modes have a Gaussian profile in the extra dimensional space. For instance, we 
depict the squared absolute values of three wavefunctions r ip T + b | 2 given by a parameter set of 
(M ab = 3 ,(ab — 0, r — i) in a direction of the torus (Rez = 0) in Fig. [U Three colored lines 
correspond to the three zero-modes labeled by I ab = 1, 2, 3. We see from this figure that the 
degenerate zero-modes are quasi-localized at points of the torus, and their peak positions are 
different from each other. These peak positions can be shifted by introducing non-vanishing 
Wilson lines ( ab ^ 0. As we said above, Yukawa couplings of the SM are expected to be given by 
an overlap integral of the zero-mode wavefunctions of a left-handed matter field, a right-handed 
matter field and a Higgs field. When their wavefunctions are localized far away from each other 
on the magnetized torus, a suppressed Yukawa coupling will be obtained. This is able to lead 
to hierarchical Yukawa couplings. Indeed, a semi-realistic pattern of the masses and mixing 
angles of the quarks and the leptons was realized without hierarchical input parameters[1]. 

Yukawa couplings can be calculated analytically by performing overlap integrals. They 
originate from the higher-dimensional gauge coupling, and the overlap integral is typically 
given by 


Y, 


IabJbcKca 9 


It 2 


dzdz (z + ( ab /M ab , r) Q J ^ (z + ( bc /M bc , r) (z + ( ca /M ca , r), 


where g is the higher-dimensional gauge coupling constant. Ref. |2J performed this integral and 
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found the result as 


Y, 


1 ab^bc^ci 


J^M ai J^M bc 
■NM ca 

X E ® 


exp 


m 


x5 


Im r 

M bc I ab —M ab J bc +M ab M bc £ 

M ab M bc M ac 

0 

lab+Jbc+Mab^iK ca ") 


> T Cab * T C be j- T Cca 

Caftlm —-h Chef™ TT-h CyjllTl — 

M ab M bc Me, 


(M bc Cab - M ab ( bc , tM ab M bc M ac ) 


( 6 ) 


for M 0ib , M bc > 0 and M ca < 0. In evaluating this result numerically, we see that the Yukawa 
couplings can be hierarchical enough to realize the ratio of masses of the top and up quarks 

(~ <?(io- 5 )). 

The largest entry of Y IabJbcKca , which is identified with a (3, 3) entry of Yukawa matrices, is 
obtained to be of 0( 1) invariably. One can easily see that, because the dominant contribution 
is given by a = 0 and l = 0 in the Jacobi-theta function (HI) . Thus, tiny neutrino mass 
{y < O(leV)) cannot be obtained with these Yukawa couplings, without heavy Majorana 
masses for the right-handed neutrinos[5j. (The top/bottom ratio can be realized by using the 
degree of freedom of tan/3 = (H u )/(Hd), where H u (Ha), represents the up- (down-) type 
Higgs held) This paper provides a way to yield the absolutely tiny Yukawa couplings in the 
magnetized toroidal compactihcations, which has a global strong suppression factor. 


3 Suppression factor in mixed SYM systems 

This section explains the basic idea for realizing a global factor to strongly suppress the Yukawa 
couplings. We consider a mixture of two SYM theories defined in different dimensions of 
spacetime, especially, we focus on a system consisting of 6D and 10D SYM theories which 
is well motivated by stable D-brane systems in type IIB superstring theory. These mixture 
configurations generically contain bifundamental representations charged under both the 6D 
and 10D SYM theories. Zero-mode wavefunctions of such bifundamentals depends on all the 
10D coordinates, but they must be localized as a point in four extra directions where the 6D 
SYM fields cannot move. 

We consider a Yukawa coupling among the two bifundamental fields and an adjoint held 
of 10D SYM theory. The two bifundamental helds should be localized at a point in the four 
direction but the 10D adjoint held has a Gaussian prohle on the magnetized tori. This Yukawa 
coupling in the 4D effective theory is determined by overlap integrals of zero-mode wavefunc¬ 
tions on each of the three 2D tori. The hrst one is given on the torus where the 6D helds live, 
and this part is represented by Eq. ([6]). The magnetic huxes on this torus will induce the hier¬ 
archical generation structure. On the other two tori, the zero-mode wavefunctions of 10D held 
is still expressed by Eq. (J3J) , but the zero-modes of the other two bifundamental representations 
are not. The point-like localization of the bifundamental helds should be expressed by a delta 
function. Introducing the delta function to express the localizations as a point seems to be also 
a sensible proposal in a D-brane picture of superstring theories)!]. Thus, the integration on the 
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C ab/Mab X 


Imz 


Figure 2: A schematic figure of the overlap integral (J7J) is shown. The blue Gaussian represents 
the wavefunction of 10D field and the red line expresses the point-like localization of 6D fields. 


torus is typically given by m 

9 J dzdz 5 T 2 (z + x )Q I ^ M - b ( z + ( ab /M ab , r) = g 0 /oi,Mo1 ’ (~x + C ab/M ab , r), (7) 

where the delta function S T 2 (z + x) must satisfy the periodicity of torus, and X indicates the 
localization point. The index I ab labels \M ab \ flavors induced by magnetic fluxes on this torus, 
and this factor is universal for generations originating from the other tori. 

This overlap (J7J) is schematically depicted in Fig. [2j This integration can take a tiny value 
when a peak position of (c) Jah ' Mab (z + ( a b/M ab , r) is far away from a point of magnetized torus 
indicated by X - That is, the distance \ X — Cab/M ab \ determines the magnitude of overlap. 
Furthermore, this suppression factor depends on the distance exponentially, and thus, we can 
obtain very strong suppression even when the value of distance \x ~ (ab/M a b\ is of 0( 1). 

We evaluate the suppression factor to verify consistency with neutrino masses, which roughly 
requires an (9(10~ 10 ) suppression[IT]. In generic (4 + 2n)-dimensional SYM theory compactified 
on magnetized tori, the net number of degenerate zero-modes of a bifundamental representation 
is given by a product of relevant magnetic fluxes on each torus, 

n 

*<*=n • 

i =1 

Three-generation structure of matter fields must originate from a single torus, and then, mag¬ 
netic fluxes on the other tori are required not to generate extra zero-modes (i.e., N ab = 3x1x1). 
Furthermore, all the generation structure of the SM must originate from the same single torus, 
otherwise the rank of Yukawa matrices is to be reduced to one and some of matter fields would 
remain massless after the electroweak symmetry breaking. 

We consider a mixture of a 6D theory compactified on a magnetized torus ( T 2 )i and a 
10D SYM theory on magnetized (T 2 ) 1 x (T 2 ) 2 x (T 2 ) 3 in the following. When all the three 
generation structure of the SM is obtained on the first torus (T 2 ), the other magnetized torus 
should not induce generations. This set-up implies M ab = 1 and I a b = 1 in the suppression 
factor (J7J) given on the second and the third tori. Such configurations can be indeed exist, and 
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we will give a concrete model with such a configuration preserving J\f = 1 SUSY in the next 
section. 

The suppression factor (JTJ) on the Uth torus is a function of two parameters, U and X 1 = 
Cab/M l ah ~ X 1 ■ Note that, r* denotes the complex structure of the i-th torus, and X * represents a 
distance between the peak position of Gaussian-like wavefunction of a 10D field and a point at 
which other bifundamental fields are strongly localized by the delta function. In the following 
analysis, we neglect the gauge coupling constant g and the normalization factor (|5]) contained 
in the function @ /,M , because their product 

gj\fi = (21m r) 1 ^ 4 e^ 2 

is always estimated as an 0(1) faetoi@ as long as the compactification scale is as high as 
Mqut = 2.0 x 10 16 . Through these, the suppression factor ((7|) is rewritten as 

f(X\r 

where we use the following property of wavefunction, 0 1,1 = 0 0,1 . 

Our numerical analysis of the function (jHJ) reveals that this factor can lead to tiny neutrino 
masses without strong fine tuning, which is shown in Table [0 The suppression can be caused on 


_ gTriXhm X i /Im r' 


E 


^7T %fir % g 2 , kHX 1 


( 8 ) 




easel 

case2 

Irn X 2 

15 

7.5 

Im t 2 

30 

15 

f 2 (X 2 , T 2 ) 

1.2 x 10~ 10 

1.5 x 10” 5 

Im X s 

0 

7.5 

Im r 3 

1.0 

15 

CO 

CO 

CO 

1.1 

1.5 x 10~ 5 

fx/ 3 

1.3 x IQ’ 10 

2.3 x 10” 1U 


Table 1: Function ()8|) is evaluated in two typical cases, which shows that the sufficient sup¬ 
pression to yield tiny neutrino masses is obtained with 0(10) values of input parameters. 


the second and the third tori, and the 4D effective couplings contain their product, / 2 x / 3 . In 
the easel, the second torus generates the tiny value of / 2 and the other is of 0(1). On the other 
side, both functions, / 2 and / 3 , contribute evenly, where the value of input parameters become 
mild. In both cases, the values of input parameters are not so sensitive. With deviations of 
(±0.5) from the sample values shown in the table, almost the same result would be obtained. 
The real parts of X 1 and r l just induce unphysical phases in this factor. 

1 We may choose a slightly large value of Imr ~ 0(10) to get a sufficiently strong suppression, but it would 

not affect on this discussion because this neglected prefactor contains only the fourth root of it. 
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4 A concrete model 


We construct a specific model of the tiny neutrino Yukawa couplings respecting the minimal 
supersymmetric standard model (MSSM), on the basis of 6D and 10D SYM theories compact- 
ified on magnetized tori. A mixture of 6D {7(1) gauge theory compactified on (T 2 )i and 10D 
{7(7) SYM theory on (T 2 )i x (T 2 ) 2 x (T 2 )3 is our starting point. We parametrize their magnetic 
fluxes as 


Mqd 


M{ 0 d 


m, 

( Mq x 1 3 

l 0 


0 

Mi x 1 2 

0 


;) 

Mr x 1 2 / 


where Mq R represents the magnetic flux on the first torus in the {7(1) gauge theory, and M[ 0D 
expresses those on the i-th torus in the {7(7) SYM theory. In Ref. [2l [TO]. conditions to preserve 
Af = 1 SUSY out of the full N = 2,3 and 4 in terms of the 4D supercharges were studied, and 
it is given by 


rn — 0, 


\ ^ M 10D = n 

^ A i 


(9) 


where A 1 is the area of the i-th torus. This form of the magnetic fluxes generically break the 
gauge symmetry as {7(7) —* U{3)c x {7(2)# x {7(2)#. The last {7(2)# gauge should be further 
broken to {7(1)#/ x {7(1)#// by introducing Wilson lines of the form 


0 0 0 \ 

Cl x 1 2 0 0 

0 0 ' 

0 0 cW 

In this parametrization, {7(7) adjoint representations are decomposed, which yield the 
MSSM fields except for the leptons. Bifundamental representations of {7(7) x 1/(1) accom¬ 
modates the charged lepton multiplets and the neutrino multiplets. Note that, the suppression 
factor is given to Yukawa couplings of the charged leptons as well as the neutrinos, but the 
distance X 1 for the each sector can be chosen independently. Thus, we can set them to give rise 
to just an 0( 1) factor for the charged leptons while realizing tiny neutrino Yukawa couplings. 

We can find a configuration of magnetic fluxes, which leads to three generation structure 
for the quarks and leptons satisfying conditions (j2j) to preserve J\f = 1 SUSY, as 


Clod — 


(Cc x 1 3 

0 
0 
0 


Mql, 

{Mb, Ml Ml) 
{M 2 c , Ml M 2 R ) 
(M 3 c , Ml M 3 r ) 


= 0 , 

= (0, 3,-3), 
= ( 0 ,- 1 , 0 ), 
= ( 0 , 0 , 1 ). 


The SUSY condition is clearly satisfied with A 1 /A 2 = A}/A A = 3. This configuration leads to 
three generations of (3, 2,1) and (3,1, 2) representation of U(3)c x {7(2)# x {7(2)#, which can be 





identified with the quark multiplets. Six generations of (1, 2, 2) representation correspond to the 
Higgs multiplets. The existence of multiple Higgs fields is a generic feature of magnetized SYM 
models, and we identify a linear combination of them with the MSSM Higgs multiplets. The 
other five combinations should be decoupled somehow. The three generations of the charged 
lepton and neutrino multiplets are correctly given by the bifundamental representation charged 
under the original 0(1) gauge group as we mentioned above. Thus, Yukawa couplings including 
the leptons have the global suppression factor (J7|). 

On this magnetized background, we can obtain 4D effective action in accordance with an 
elegant way[9j, and there, specific spectra of the quarks, the leptons, and the SUSY particles 
can also be calculated. Although we will not do that because that is not our purpose in this 
paper, we are able to easily infer that similar spectra to those obtained in Ref. [!.] can be 
realized except for the neutrino sector. 


5 Conclusions and Discussions 

We have studied strong suppression due to quasi-localization of wavefunctions and realized 
tiny neutrino Yukawa couplings in magnetized toroidal compactihcations of SYM theories. In 
a system of two SYM theories defined in different dimensions, such as, 6D-10D SYM theories, 
bifundamental representations of the two gauge groups must be localized as a point in extra 
four directions, which should be represented by the well-defined delta function. As a result, 
their 4D effective coupling with a 10D held contains a factor induced by overlap integrals of 
the delta function and a Gaussian wavefunction of the 10D held in the four directions. When 
a point on which the delta function will not vanish is equal to the peak position of Gaussian, 
the integral induces an 0(1) factor. However, for their nonvanishing distance, the value of the 
factor is reduced exponentially. Our numerical analysis shows that the 0(10) of distance leads 
to an 0 ( 1 ( 0 10 ) of suppression, which is strong enough to yield tiny neutrino masses consistent 
with experimental and observational data. 

A specihc SYM system and its configuration have been discussed in Section 4. We consider 
a mixture of a 6 D 17(1) gauge theory compactihed on magnetized (T 2 )i and a 10D 0( 7) SYM 
theory on ( T 2 )i x (T 2 ) 2 x (T 2 ) 3 . The 0(7) gauge group is broken in part by magnetic huxes 
to U(3)c x 0( 2)l x U{ 2) r . The decomposed 0( 7) adjoint held accommodates the MSSM 
contents except for the lepton multiplets. The lepton multiplets assigned into the bifundamental 
representation charged under the original ( 6 D) 0(1) gauge group. Thus, Yukawa couplings for 
the charged leptons and the neutrinos are determined by an overlap integral of three Gaussian 
wavefunctions on the hrst torus (T 2 )i, and by overlap integrals of the delta function and the 
Gaussian on (T 2 ) 2 and (T 2 ) 3 . Since the distance X 1 between the peak of the delta function 
and the Gaussian can be controlled independently for the charged leptons and the neutrinos, 
neutrino Yukawa couplings can be strongly suppressed while just an 0(1) factor is given rise 
to for the charged lepton sector. 

This very strong suppression could be applied to other couplings to forbid dangerous pro¬ 
cesses, such as, the proton decay and the flavor violation. In magnetized toroidal compactihca¬ 
tions, higher-order couplings can also be calculated analytically[ 12 ], and this strong suppression 
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mechanism would also work in such higher-order couplings as well as Yukawa couplings. Thus, 
we expect that a wide variety of its applications are possible. We will construct models focusing 
on other phenomenologies and analyze them elsewhere. 
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